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Abstract 

We postulate a master equation, written in the language of creation and annihilation operators, as 
a candidate for unambiguous quantum mechanical description of unstable particles. We have found 
Kraus representation for the evolution driven by this master equation and study its properties. 
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I. INTRODUCTION 


One of important difficulties in quantum mechanical description of unstable particles is 
irreversibility of time-evolution. The complete system consists of decaying particle as well 
as of decay products. Only this complete system undergoes unitary evolution, described 
by quantum field theory. However, in many applications of quantum mechanics (e.g., in 
analysis of correlations experiments) we would like to neglect the evolution of decay products 
and consider the decaying particles only. This is usually achieved by introducing a non- 
Hermitian Hamiltonian, as it was done in the classical works of Weisskopf and Wigner 
00- The non-Hermitian Hamiltonian, however, leads inevitably to the non-conservation 
of trace of the density operator of the system. Although such description gives decreasing 
probability of detecting the particle, it does not provide unambiguous way of calculating the 
probability of finding the system consisting of a few such particles in a given state after the 
measurement (see e.g. [3| for a discussion of other ambiguities caused by various approaches 
to the description of such a system). Indeed, one needs to use probability theory rather 
than quantum mechanics for this purpose. Since the calculations of conditional probabilities 
are extremely important for analysis of correlation experiments, especially of those done for 
systems of neutral kaons or H-mesons [ThMTsI] . it would be desirable to find a quantum 

mechanical description of decaying particles preserving unit trace and positivity 16|-22| of 

], 24j]) show that there is still a 


the density operator for the system. Recent papers (e.g. 


great interest in the unambiguous quantum-mechanical description of neutral kaon system. 
This is the point where the theory of dynamical semigroups and open quantum systems 

[25H27 


can be helpful. Let us recall (cf. [28H30 


that the dynamical semigroup in the 
Schrodinger picture is a one-parameter family of linear maps A*, acting on the space of 
trace class operators on Hilbert space of the system, preserving for every t > 0: (i) positivity, 
(ii) trace, (iii) strong continuity and such that (iv) A£ A* = A * 1+t2 for every t\. t 2 > 0. These 
properties can be translated into Heisenberg picture as requirements for the map A t acting 
on the space of bounded operators on the Hilbert space of the system, which for every t > 0: 
(i) preserves the positive cone, (ii) leaves the identity operator invariant, (iii) is continuous 
on states in the trace-norm sense, (iv) is normal and (v) A tl A t2 = A tl+t2 for R, f 2 > 0. 

The idea that the theory of open quantu m sy stem would be useful for the description 
of unstable particles appeared quite early [3ll^33| (see also 34] for a review). Recently, the 
open quantum system approach was also applied to the systems of particles with flavour 


oscillations (like in the case of neutral kaons) [35], [36|, and it has been used successfully in 


the description of EPR correlations and evolution of entanglement in K°K° system [sul . 38 

Here, we follow the approach presented in [35, 37, 38jj. However, in these works the 
considerations were restricted to systems of at most two particles, and transition from one- 
particle to two-particle theory was done by means of tensor product construction. In the 
present paper we will show that it is possible to describe systems with arbitrary number of 
particles using the second quantization formalism, which is the most natural language for 
system with varying number of particles. Moreover, such approach would be an advantage 
if we study the behaviour of the system from uniformly moving or accelerated frame, due to 
the well established transformation properties of annihilation and creation operators. 

The paper is organized as follows. In Sect. HQ we postulate a master equation in 
Schrddinger picture for a single kind of free particles, and then find the solution of this 
equation in the form of the Kraus representation of the evolution of the density operator 
of the system. The next section is devoted to the Heisenberg picture of the evolution of 
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the same system. In Sect. El we analyse the system of particles of different types and the 
flavour oscillation phenomenon. 


II. SCHRODINGER PICTURE 


it was shown that the time evolution of a free unstable scalar particle can be described 


In 

by a master equation in the Lindblad-Gorini-Kossakowski-Sudarshan form 
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^ + Lp(t)L *, (1) 

where 

H = m\ 1)(1|, L = V f|0)(l|, K=~L^L. (2) 

Here |1) denotes the state of presence of the particle and |0) denotes the state of its absence; 
m is the mass of the particle and Y is its decay width. Despite the fact that the state |0) 
is usually called vacuum, it is not the vacuum in the sense of quantum field theory, but 
rather in the sense used in [39], i.e., it is an absence of a particle. This equation leads 
to the probability density of finding the particle evolving according to the Geiger-Nutall 
exponential law. 

However, the most natural quantum-mechanical description of systems with variable num¬ 
ber of particles is the second quantization formalism. For systems governed by (JTJ) the tran¬ 
sition to second quantization is straightforward, since the operators ([2D can be interpreted 
as the vacuum-one-particle sector of the second quantized operators 


H = ma l d = mN, L = VTa, (3) 

where a and a* are bosonic annihilation and creation operators, respectively: 

[a, a f ] = 1, a|0) = 0, \n) = -j=- 10) ; (4) 

Vffi 

vectors | n) form the so-called occupation number basis. With this substitution we have 

K = --YN. 

2 

If we substitute (J2D into (dD, we arrive to the master equation in the form 



£*(A», 


(5a) 


where 


£*(p) = - im[N , p\ - —{2V, p} + rapa f . 

In the above we can recognize the equation introduced and studied in 
equations lead to evolution given by quasi-free semigroups, see eg. |3i 
studied in the context of quantum optics 


148-51 


9,31, 

32 

,34] 

o 

- 

and 


(5b) 


Similar 
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If we write down explicitly annihilation and creation operators in occupation number 
basis, namely a ki = y/k + 1 4+u, + I4,z+i, so iV fcZ = k5 ki (k, 1 = 0,1,...), then we 

can view (ED as the following infinite system of equations for matrix elements of the density 
operator A*p = Jf kl Pki(t ) |fc) (Z|: 


dpkiif)_ 

dt 


i(k — l)m — \(k + l)T p k i(t) + 


\!(k -\-!)(/ + l)I> fc+u+ i(t), 


( 6 ) 


for k, l = 0,1,.... Thus we get infinite, in principle, system of linear differential equations 
of first order. 

Notice, that the system (ED seems to be highly non-trivial — the solution for p k i(t) depends 
on a solution for pfc+i,z+i(£), what apparently leads to infinite chain of dependencies. What 
makes the system (ED solvable is the proper choice of initial conditions. Indeed, for every 
reasonable initial physical state, the number of particles must be finite, so all matrix elements 
of p( 0) corresponding to higher number of particles must vanish. Mathematically, it means 
that there exist indices r and s, such that 


Pki (0) = 0 for k > r and l > s . 


(7) 


One can easily check that the system (ED with initial condition (ED gives us the well posed 
Cauchy’s problem. 

Now, instead directly solving the equation (ED for some interes ting choices of initial state 
we concentrate on finding and studying the Kraus representation 52] of the evolution of the 
system. 

Although Kraus operators for the evolution of the density operator governed by the 
master equation (ED was found in 53|, here we give its another formulation, now written in 


terms of annihilation/creation operators. It can be easily checked that these two choices of 
Kraus operators coincide up to the phase factors. Despite this, we give the formal proof that 
proposed Kraus operators lead to the evolution of the system undergoing the equation (ED, 
since we will employ the technique used in the proof later on. 

Proposition 1. If for k = 0,1,.. 


E k {t) = 


—iMt 


yfk\ 


e~ rt a 


where M — (m — |T)7V, then 


A ; P = '£E k (t) P El(t) 


( 8 ) 


(9) 


k =0 


is the solution of the master equation (ED , where p is the density operator given at initial 
time t = 0. 


Proof. First, let us note that for any k 

E k (t)d = e { - im+ ^ )t aE k (t), (10) 

what follows immediately from canonical commutation relations. Using (fTUD we can show 
by straightforward calculation that 

dEo(t) 


dt 

dE k (t) 

dt 


= ~iME 0 (t ), 
= — iME k (t) + 


v /^ Fe (im+ ±r)t 
2\/l - e~ rt 


aE k _ x (t) 


(11a) 

(lib) 
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(for k = 1,2,...). 

Next, one can easily check that the following recurrence relations hold for k — 1,2,... 


E k (t) = Vl e r V im +i r > t a E k ^(t) 


( 12 ) 


Combining (fT2|) and (lllbl) we can write (HID in the form 


dE k (t) ( - k Te~ rt \ . , 

~dt~~ \~ lM + 2T~^) Ek{t) ' 

Now, let us compute the time derivative of the density operator A*p given by 


(13) 


dt , 


.dEtm 


k =0 




p„—rt oo 

M(Kp) - (AMM»] + —— FJ E kE t (t)pEl(t) . (14) 

1 6 k =1 


Taking into account the last term in (UTl) can be written as 

rt. oo oo 

, _ _ rt J2 kE k (t)pEl(t) = TaJ2 = ra(A t »a f . (15) 

1 e fe=l k=1 

Thus, the density operator obeys the master equation (EJ. 

To complete the proof, we have to show that at the time t — 0 the density operator A *p 
given by (P) is p. This is trivial point, because obviously _E 0 (0) = id and E k ( 0) = 0 for 
fc = l,2,... ^ □ 

It is easy to see that after writing out annihilation operators in occupation number basis 
the Kraus operators ([SD differ only by phase factors from those found in pki 54], These 
phase factors become important when you try to study a flavour oscillation phenomenon 
(see Sect. ED. 

Proposition 2. If E k {t) are given by (I8JL then 


E E l(t)E k (t) = id . 


(16) 


k =0 


Proof. We start with the observation that for any element of occupation number basis | n) 
and any non-negative integer k 


n\ 


a fe \n) = { V (n-k)! 

0, 


\n — k) , n > k , 
n < k , 


(17a) 


and 


(at) fc |n) = |n + k) . 


n\ 


(17b) 
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Since 


4(t)E k (t) = i (1 - e' r, )‘ (at) t e-™a‘, 

then, from (TT7T) . for any element of the occupation number basis 

El(t)E k (t) | n) = Q (l - e ~ vt ) k e- {n ~ k:)Vt \n) , (18) 

when n > k, and E k (t)E k (t ) | n) = 0, when n < k. Thus, 

oo n n / \ , 

E 4(‘)^(f) |n) = E El(t)E„(t) |n) = E J (l - I") = I") ■ (19) 

k =0 k =0 A:=0 \ K J 

Since YfkLo £* (£)-£*(£) acts as identity on any element of the basis, it must be the identity 
operator. □ 

Proposition 3. Vacuum state |0) (0| is stable under the evolution given by (jHJ) and (ED- 
Moreover, lim t _ ) . 00 A*p = |0) (0| for any density operator p. 

Proof. Indeed, E 0 (t) |0) = |0) and E k (t) |0) = 0 for k = 1, 2,..., so the density operator 
|0) (0| is stable during the time evolution. 

For the proof of the second statement, one can find that 


E k (t) | n) 



when n > k, and vanishes otherwise. Next, observe that for T > 0 


( 20 ) 


jto e -W-‘X (^^)‘= j 

Loo, 

n > k , 
n = k , 
n < k , 

(21) 

SO, 




, o, n V k , 

Inn E k (t) ra) = < 

10), n = k , 


(22) 

and, consequently, 




lim E k (t) |n> <n'| E[(t) = |Q) (Q| _ 

n n !, 
n — n'. 

(23) 

Thus, for any density operator p, lim^oo A*p = tr(p) 0) (0 = 

= |0) (o|. 

□ 


Explicit solutions of ([TD with operators (J3J) can be deduced from the relation 


minlnjn'} 


A * \n) (n'\ = Y 


n\ in 


=o \\ k J\ k , 


0 —im(n—n r )t 


X e -\T{n+n'-2k)t ^ _ e ~Tt^ k ^ ^ _ fc | _ ( 2 4) 
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If we impose the superselection rule which forbids the superpositions of states with different 
number of particles, then the density operator for a system consisting of at most n particles 
is of the form 

n n 

Kp=Y,Pk{t)\k) {k\ , J2Pk(t) = 1- 

k =0 k =0 

Thus, it is enough to solve the equation (P) for an initial state of the form p = \n) (n\ for 
n being some non-negative integer (arbitrary, but finite), because any density operator for 
the initial state is a linear combination of such states. 

If the system is initially in the n-particle pure state, p = | n) (n\, then the solution of the 
equation (JTj) is 

A*p = E (l) e ~ {n ~ k)rt (! - e~ rt ) k \n - k) (n - k\ . 
k =0 \kj 

The average number of particles changes in time as 


(N(t)) = tr (A ;p)N 


= ne 


-Tt 


We have thus recovered the Geiger-Nutall exponential decay law. It is worth noting that 
the probability that at a time t one finds exactly k particles from initally n ones, has a 
binomial distribution B(n, e~ rt ) with probability e~ vt of finding a single particle, as it can 
be expected. 


III. HEISENBERG PICTURE 

In Heisenberg picture, master equation for the evolution of an observable is of the form 

4A= C{A t n ), (25a) 


where 


£(fi) 




(25b) 


Note that 

^{[L\ Q]L + L f [n, L\} = {K, G} + LHiL , 

but the form used in (I25b|) is usually more convenient when performing calculations in the 
Heisenberg picture involving creation and annihilation operators. 

Having a family of Kraus operators ([HD, the evolution of observable Q can be written as 
the series 

OO 

A t n = Y,El{t)nE k {t). ( 26 ) 

k =0 

This representation is especially useful if we can find the decomposition of the observable 
into its matrix elements in occupation basis: 


fl = E^n.n' \n) (n' 

n,n' 


(27) 
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A t | n) (ri 


E \ 

k =0 \ 



( n ' + k\ im(n—n')t 

{ n> ) 6 


x e -\nn+n')t ^ _ g-r t y | n + ^ + fc | _ (28) 

Using projectors onto n-particle states, n n = |n) (n|, the last equation can be rewritten in 
a more convenient form 


Atfl n 



E 




(29) 


Proposition 4. Hindoo A t n 0 = id. 
Proof. From (f2Ul) it follows that 


°° h 

A4o = £ (i - e- R ) n fc> 

fc =0 

oo 

so lim^oo A t n 0 = £ n fc = id. □ 

k =0 

Physically, Proposition [I] tells us that after substantially long (mathematically infinite) pe¬ 
riod of time, the probability of finding vacuum reaches one, irrespectively of the state of the 
system. In other words, at infinite time all the Fock spaces collapse to the vacuum subspace. 

The evolution of creation and annihilation operators can be easily find with use of rela¬ 
tion (HUD : 


A t a = e _( * m +l r )*Q } (30a) 

A t a f = . (30b) 

Moreover, it is easy to check that in this case A t N = A t af A t a (what, in general, does not 
hold). Indeed, using (fIU)l we have 

OO / OO \ 

A t N = £ El(t)a)aE k (t) = £ E^E k (t) = A t a^A t a . (31) 

k= 0 \fc=0 ) 

Consequently, the evolution of the particle number observable is 

A t N = e~ rt N- (32) 


we can get this result by solving (I25aj) for N, too. 

It is easy to find the mean number of particles for a given state with help of (1321) . Here, 
we consider two examples: the pure state of exactly n particles and a coherent state with 
given mean number of particles n. 

Example 1. If the system is in the pure state of n particles, then the mean number of particles 
is simply 


<JV(()) = ne- n . 


(33) 






Thus we get the exponential decay law again. Time evolution of the probability of finding 
exactly k particles follows from (I2U1) and reads 


Pn(k,t) = 



1 


i.e., it is given by the binomial distribution B(n,e rt ). 

Example 2. Let us assume that the system is in a coherent state |a), 


a\a) — a\a) , 


(34) 


a E C, i.e. 


then 


a) 



OO 


E 



\k) , 


(N(t)) = ne ~ Tt , 


(35) 

(36) 


where n = |cr| 2 is the mean number of particles in the coherent state |a). Probability of 
finding exactly k particles evolves in time according to 


Pn(k,t) = i (ne r ') t e “ " , 


(37) 


which is the Poisson distribution P(ne~ Tt ). 

Let us note that, if we consider the state being a mixture of fc-particle states with prob¬ 
ability that fc-particle state occurs given by the Poisson distribution with mean number of 
particles n, i.e., 


p = E 

k=0 


oo e n n k 


k\ 


I k) (k | , 


(38) 


then the mean number of particles in this state and probability of finding exactly k particles 
are given by the formulae and (TTT7|i . respectively (despite the fact, that in this case we 
must find the traces of the product of observables with the density operator). 


IV. PARTICLES OF DIFFERENT TYPES 


Let us consider a system of particles of r different types (or carrying a quantum number 
with r possible values), each type with mass rrij and width Tj for j = 1,,r. For such a 
system we have 

[Sj, a,k]zf = 0 , [a,j, a k ]if = 5jk , (39) 

for j, k — 1 ,... ,r, where [•, denotes commutator/anti-commutator, respectively, and anti¬ 
commutators apply only if both j th - and £; th -particles are fermions. The states spanning the 
occupation number representation are generated from the vacuum state via the formula 


^11 ^2 j • ■ ■ j ^ r ) 


(ai>-(air-■■(&;)"' r 

V'n 1 !n 2 ! • • • n r \ 


(40) 


where we identify |0, 0 ,..., 0) = |0). 
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The master equation for the system takes the following forms 


Up = + {it, A,» + £ 4(VV)4 , (41a) 

ar j =1 

^A t n = i[#,A t n] + X!{[Lt j A t n]L i + Lt[A t n,4]}> ( 41b ) 


in the Schrodinger and Heisenberg picture, respectively, where H is the Hamiltonian of the 
system and 


4 = V / 4'V K = - 1 -jZL]L i , M — H + iK. 

Z j =1 


(42) 


If [M, hj] = —(jrij — |T j)aj for j = 1,..., r, then we can easily construct the Kraus operators 
solving QHD 


E k {t) = e 


—iMt 


n 


f V 1 — e 1 j!e> : ) 


(43) 


k\ 

k\ H - ” "~\~kj' := k 


where the product is taken over all possible partitions of k into exactly r addends, such that 
k\ + A )2 + • • • + k r — k, where 


kj G N 0 , 
kj G {0,1} , 


j th -particles are bosons , (44a) 

j th -particle is a fermion . (44b) 


for j = 1,... ,r. 

Example 3. Let us consider the evolution of a system of two flavour particles (e.g., particles 
and their anti-particles). We denote the creation operators for these particles by a\ and a)>. 
The basis for the system is built up from the states of the form 


ni,n 2 ) 


(a t 1 ) ni (a^) n2 

Vni\n 2 \ 


(45) 


Let these states be the common eigenstates of two observables, N = a\ai + a 2 a 2 (number of 
particles) and S = a\ai — a 2 a 2 (say strangeness or lepton number), i.e., 


N\m, n 2 ) = (rai + n 2 ) |ni, n 2 ) , 
S\ni,n 2 ) = (ni - n 2 ) |ni, n 2 ) . 


If the states (H51i are not eigenstates of the time evolution the phenomenon known as the 
flavour oscillation may occur. 

To describe such a situation, let us assume that the Hamiltonian and Lindblad operators 
for the system are of the form 


H 

L\ 

L 2 


rriic\ci + m 2 c 2 c 2 , 

\/r\ci, 

^2, 


(46a) 

(46b) 

(46c) 
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where cj,c 2 are connected with aj,a 2 by unitary transformation: 

at _ Jx ( Priq I -i(</>-y )/2 


t *t 


c = e A e 


cos -a[ + e 


sm -a 2 j , 


at _ Jx I _J(<t>-4>)/2 Q ; n , _-i(0-H/O/2 


Co = e 


sin -a[ + e 


cos - a, 
2 2 . 


(47a) 

(47b) 


Since M = (Vp — fTi) cjci + (m 2 — |T 2 ) c 2 c 2 , we can easily find the evolution of cf 


hcj 

Using (147)1 we get the evolution of a 3 : 


A,c, = e ~( im ‘ + i r ‘)% , j = 1,2. 


(48) 


AfSi = -e 


— imi + 


5 r l> 


ai(l + cos 9) + a 2 e^ sin 9 


+ -e ( ,m 2 + 2 r 2 )f ai(l — cos6>) — a 2 e^sind 
a 2 (l + cos 9) — aie ~sin 9 


A t a 2 = i e -( imi+ 5 ri )' 


+-e ( im2+ 2 r 2) 4 a 2 (l — cos6*) + hie ^sind 


(49a) 


(49b) 


The time evolution of the observables can be obtained either by solving (141bl) or directly 
from relations fl49l) . using argumentation analogous to (13111 . For example, for the number of 
particles we get 


A t N = 


e Flt + e T2t 


N + 


e -rn _ e -r 2 t 


[S cos 6 + Q + sin 9] 


2 2 

where Q + = + a 2 aie - *^, so the mean value in the state |ni, n 2 ) is the following 


(50) 


(N(t)) = 


3 —ri t _|_ g—r 2 t 


■{ni + n 2 ) + 


g—rn e -r 2 t 


2 v 2 

and is depicted in the Fig. [0 

Similarly, for the strangeness (or lepton number) we get 


(ni — n 2 ) cos 9 


(51) 


A ,S = 


g—Tii g—r 2 t 


N cos 9 + e rt sm(Amt)Q- sin 9 


+ 


e Fl4 + e V2t 2 


cos 2 9 + e r * cos(A mt) sin 2 9 


S 


+ 


g Tit 4- g r 2i 

---e cos(Amt) 


Q+ sin 9 cos 9 , (52) 


where Q _ = i(ala 2 e^ — a 2 aie l?i ), T = |(Ti + T 2 ) and Am = m 2 — mi. The mean value of 
this observable in the state |ni,n 2 ) is 


(Sit)) = 


g—Tit_g—r 2 t 


(ni + n 2 ) cos0 

+ 


e rit _i_ e r 2 t ^ 


cos" d + e 14 cos(Amt) sin" 9 


(n i - n 2 ) (53) 
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FIG. 1. Number of particles for system with two flavours for mixing angles 6 - 
right to left) with n\ = 2 and n 2 = 1, and ri < T 2 (time unit is t = 1/r). 

= 0, (from 

and is shown in the Fig. [21 

Let us pay our attention on the two extreme cases: 9 = 0 (no flavour mixing) and 9 — | 
(maximal mixing). For 9 = 0 we have q = q, so the time evolution of the observables is 

A tN = \e~ v '\N + S) + - S) , 

(54a) 

A t S = \e~ rit (S + N) + ^e~ r2t (S - N) . 

(54b) 

Their mean values in the state \ni,ri 2 ) are 

(N(t)) = e~ Vlt n i + e _r2t n 2 , 

(S(t)) = e~ llt ni — e~ r2t n 2 . 

For 9 — <j) — 27t, ip — 7T and x — we have 

(55a) 

(55b) 

c\ = ^|(a f i + a\) , 
c 2 = ^(«i - a\) , 

and the time evolution of the observables is given by 


A t N = i (e"™ + e" r2 *) iV + 1 (V^ - e“ r2t ) Q + , 

(56a) 

A t /S = e _r * cos(A mt)S + e _It sin(Amf)Q- • 

(56b) 
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0 1 2 3 4 5 


t 

FIG. 2. Strangeness of system with two flavours for different values of mixing angle with n\ = 2 
and n 2 = 1, and Fj < T 2 (time unit is r = 1/r). 


The mean values of these observables in the state |ni,ri 2 ) are 

(N{t)) = ^ (e~ Tlt + e _r2t ) (m + n 2 ), (57a) 

( S(t )) = e~ rt cos(Amt)(ni — n 2 ), (57b) 

so we get oscillations of the quantum number S. It is worth noting, that for the particles 
such as K or B mesons we can use this result only as a first approximation, since for these 
particles the transformation which “diagonalizes” the master equation is non-unitary due to 
UP-violation. For the sake of brevity, we do not discuss the violated CP-symmetry here 
but preliminary calculations show the agreement with values for masses and life-times of 
neutral K or B mesons estimated on the basis of traditional Wigner-Weisskopf approach. 


V. CONCLUSIONS 

We have analyzed a class of master equations built up from creation and annihilation op¬ 
erators which generate dynamical semigroups that can describe the exponential decay and 
flavour oscillations for system of many particles. We have shown, in this case, how this 
dynamical semigroup can be written in the Schrodinger as well as Heisenberg picture. This 
allowed us to choose the picture which seems to be more convenient for the description of 
the system under consideration. Moreover, we have found the solution for a free particle 
master equation in the form of Kraus representation in the language of annihilation and 
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creation operators. Although, this Kraus representation is given by an infinite series, in the 
Schrodinger picture it reduces to a finite sum, whenever the initial state has a finite number 
of particles. On the other hand, in the Heisenberg picture the commutation relations be¬ 
tween observables and Kraus operators sometimes allows us to find the observable evolution 
in closed form without explicit summation of the series. 

Notice that if we cut the presented approach to the one-zero particle sector we get the 
theory given in 35] (neglecting the decoherence). 

In the present paper we restrict our analysis only to states labeled by a discrete index, 
and not by continuous parameter (like e.g. momentum). Despite the fact that introducing a 
continuous parameter causes creation and annihilation operators to become operator-valued 
distributions, it seems to us that the approach introduced here should also be applicable. 

We left open the question whether it is possible to apply our approach to describe the 
processes other than exponential decay, like e.g. decoherence or different decay laws. The 
preliminary investigations suggest that there exists a positive answer. 
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